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We study topological defects as inhomogeneous (localized) condensates of particles in Quantum 
Field Theory. In the framework of the Closed-Time-Path formalism, we consider explicitly a (1 + 1) 
dimensional Xtp 4 model and construct the Heisenberg picture field operator ip in the presence of 
kinks. We show how the classical kink solutions emerge from the vacuum expectation value of such 
an operator in the Born approximation and/or A — *• limit. The presented method is general in 
the sense that applies also to the case of finite temperature and to non-equilibrium; it also allows 
j^JT)' for the determination of Green's functions in the presence of topological defects. We discuss the 

[""t ' classical kink solutions at T 7^ in the high temperature limit. We conclude with some specu- 

lations on the possible relevance of our method for the description of the defect formation during 
symmetry-breaking phase transitions. 
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I. INTRODUCTION 



In recent years much attention has been devoted to the study of topological defects and to the issue of their formation 
in the course of symmetry-breaking phase transitions (l],^]. This interest arises essentially for two reasons. On one 
hand, topological defects are naturally inherent to many dynamical systems ranging from cosmology to condensed 
matter physics, and thus they provide a bridge between processes which manifest themselves at very different energies 
Q_i! and time scales j^H. On the other hand, the study of the formation of defects during phase transitions offers an 
important way to understand better the underlying non-equilibrium dynamics of quantum fields H]?]] . 

The basic phenomenological picture of defect formation is presently understood via the Kibble-Zurek mechanism 
jl],Q: In many systems, when the phase transition associated with a spontaneous symmetry breaking takes place, 
some regions of space may remain trapped in the initial symmetric phase. These regions are called topological defects. 
In this way cosmic strings may have been formed in the early universe; vortices in superfluids and superconductors, 
and textures in some nematic liquid crystals are also created in a similar fashion. A full understanding of the process 
of defect formation in the context of Quantum Field Theory (QFT) is however still far to be reached, although recently 
there has been much progress in this direction [|fj|-|l2|. There are two fundamental points which make difficult such a 
task: the emergence of the macroscopic (collective) behavior associated with the defects in question resulting out of 
the underlying microscopic (quantum) dynamics, and the intrinsic off-equilibrium character of the processes - realistic 
phase transitions - in which the defects are created. 

In the present paper, we address the first question, i.e. we show how a self-consistent description of topological defects 
can be achieved in QFT in a way that is suitable for the study of their formation in non-equilibrium processes - task 
which we plan to tackle in our future work. We are inspired by an approach to the description of extended objects 
(solitons) as condensates of particles in QFT which was developed in 70's by Umezawa, Semenoff et al. |b]|l4|]. In their 
approach, known as the boson transformation method, the usual order of argumentation is reversed: solitons emerge 
naturally as macroscopic objects with the full quantum origin. The classical soliton solutions are then recovered in 
the H — ► limit. Here we revisit this method and reformulate it in the framework of the Closed-Time-Path (CTP) 
formalism [ p"5| . There are several advantages in doing this. First of all, we obtain a closed functional form for the 
Heisenberg picture field operator in the presence of defects. This is to be compared with the original approach in 
which such an expression is obtained recursively via the Yang-Feldman equation. The most important feature of the 
presented formulation relies however on the fact that the CTP formalism allows for a unified treatment of both the 
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zero temperature and the finite temperature cases. It is also suitable for the treatment of non-equilibrium situations: 
in this case, as it is well known, the use of CTP is essential jl6). Last but not least, our formulation allows to 
determine Green's functions in the presence of defectSj_which can be useful for example for setting up self-consistently 
the Cauchy data for a non-equilibrium evolution Jn],|l8| . 

In this paper, we consider a toy model system, namely the (1+1) dimensional Aip 4 theory, both at zero and finite 
temperature. We construct the Heisenberg field operator in the kink sector and recover the known (classical) kink 
solutions in the classical limit as vacuum expectation values of the field operator. We do this in order to set the 
mathematical framework and explain techniques which we are going to "loosely replicate" in our further work, where 
more realistic systems will be considered |H| . 

An interesting consequence of the present method are the recurrence relations we arrive at in Section III, which reveal 
an intimate connection between the (non-exact) solvability of the (1 + 1) Aip 4 theory and the solutions of Cauchy- 
Marley's functional equations. In this paper we do not consider quantum corrections to the classical soliton solutions, 
although the method here presented allows for a systematic study of this aspect (see comments in Section III) , which 
will be treated elsewhere. 

We then go one step further and show how to extend the treatment for a system in thermal equilibrium saturated 
at some temperature T. We show that in the high T limit it is indeed possible to get analytical solutions for the 
(inhomogeneous) order parameter. A delicate limit h — > at finite T is discussed and compared with other existing 
treatments [^(| . The results obtained are of general validity and can be applied, for instance, to more realistic higher 
dimensional cases for the study of restoration of symmetry. 

The plan of the paper is as follow: In Section II we briefly review the boson transformation method and the Haag 
expansion which are the basic ingredients of our method. As we aim to study topological defects in both equilibrium 
and non-equilibrium media we formulate the whole approach in the CTP formalism. 

In Section III we consider the zero-temperature (1+1) A?/> 4 theory in the Goldstone phase. We construct the classical 
kink solutions from the vacuum expectation value of the Heisenberg picture field operator, in the Born approximation 
and/or the A — > limit. These limits are realized via residue calculations of certain functional relations. 

In Section IV we extend our study to finite temperature. To outline the basic strategy we consider the case of large 
T: this allows us to obtain in a closed form the solution for the (thermal) order parameter. The high temperature 
analysis reveals that there cannot be any phase transition at high T but, instead, the system is plagued by spontaneous 
fluctuations exemplified by the non-analytic behavior in the order parameter. We also show explicitly that one cannot 
talk about the loopwise expansion as being an expansion in h at finite temperature. This essentially distinguishes our 
approach from the ones where tree approximations at finite T are considered as equivalent to classical limit. 

We conclude in Section V with some comments and speculations on the possible relevance of the presented method 
for the description of defect formation during phase transition. We also discuss the role of the shift function (used for 
a construction of the inhomogeneous coherent states) as a possible alternative tool for the classification of topological 
defects. 

In the Appendices are clarified some finer mathematical manipulations needed in the main body of the paper. 



II. DESCRIPTION OF SOLITONS AS INHOMOGENEOUS PARTICLE CONDENSATES 

It has been recognized since long time that due to the breakdown of von Neumann's uniqueness theorem |2l| l the 
structure of Quantum Field Theory is extremely rich, allowing for a host of different (mutually unitarily inequivalent) 
Hilbcrt spaces for a given dynamics of Heisenberg fields. The choice of the Hilbert space is equivalent to fixing the 
boundary conditions for the (operatorial) Heisenberg equations. It should be born in mind that a particular choice of 
boundary conditions specifies the observational particle (or quasi-particle) content of the theory. In the perturbative 
Lehmann-Symanzik-Zimmermann (LSZ) approach, the boundary conditions are injected via the asymptotic fields 
(most commonly via the in-fields), and the functional relation between them and the Heisenberg fields is then known 
as Haag 's map. One practically realizes such a map expanding the Heisenberg fields into a (functional) power series 
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of the asymptotic fields - the Haag 's expansior^ 

It is worth to stress that Haag's expansion is not an operatorial relation, but holds true only in a weak sense, i.e. 
when matrix elements are considered. As we shall see in Section II. B, this is due to the fact that the Hilbert spaces 
for the Heisenberg fields and asymptotic fields are mutually unitarily inequivalcnt. 

A. The "boson transformation" method 

Starting from these basic observations, an approach to the description of extended objects (solitons) as condensates of 
particles in QFT was developed in the seventies by Umezawa et al. fb^Jl4[] • In this approach solitons emerge naturally 
as macroscopic objects (inhomogeneous condensate) with the full quantum origin, provided one chooses the "right" 
Hilbert space for the asymptotic fields. The classical soliton solutions are then recovered in the Born approximation. 

Let us here briefly recall the main lines of the boson transformation method. To avoid unnecessary difficulties we will 
ignore, temporary, the renormalization problem. Let us consider the simple case of a dynamics involving one scalar 
field tp satisfying the equation of motion: 

A(d)ip(x) = J[iIj](x), (1) 

where A(d) is a differential operator, x = (t,x) and J is some functional of the ip field, describing the interaction. 
Let us now denote the asymptotic field by (p(x), satisfying the equation 

A(d)<p(x) =0. (2) 

Equation ([!]) can be formally recast in the following integral form (Yang-Feldman equation): 

i/i{x) = ip(x) + A-\d) * J[ip]{x) , (3) 

where * denotes convolution. The symbol A -1 (<9) denotes formally the Green function for the (p(x) field. The precise 
form of Green's function is specified by the boundary conditions. Eq.(^|) can be solved by iteration, thus giving an 
expression for the Heisenberg field ^>{x) in terms of powers of the f(x) field; this is the Haag expansion (or "dynamical 
map" in the language of refs. ]l3| , p^ |), which might be formally written as 

= F [x; ip] . (4) 

As already remarked, such an expression is valid only in a weak sense, i.e. for the matrix elements only. This implies 
that Eq.(||) is not unique, since different sets of asymptotic fields (and the corresponding Hilbert spaces) can be used 
in the construction. Let us indeed consider a c-number function f(x), satisfying the same free equations of motion: 

A(d)f(x) = 0, (5) 

then the corresponding Yang-Feldman equation takes the form 

^{x) = <p + / + A- x (9) * J[ip f ]{x) . (6) 

The latter gives rise to a different Haag expansion for a field ipf (x) still satisfying the Heisenberg equation (Q) : 

lS(x) = F[x;<p + f\. (7) 

The difference between the two solutions ip and ip* is only in the boundary conditions. An important point is that 
the expansion Eq.(Q) is obtained from that in Eq.(^J), by the space-time dependent translation 

<p{x) -» <p(x) + f{x) . (8) 



1 The Haag expansion is often denned as a power series of a normal ordered product of the asymptotic fields The 
connection between the definition we use and the one in |E3| can be established via the operatorial Wick's theorem. 
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Eqs.fH) and (g) express the essence of the so called, boson transformation theorem J13[ : the dynamics embodied in 
Eq.(0l, contains an internal freedom, represented by the possible choices of the function f(x), satisfying the free field 
equation (JsJ) . Also observe that the transformation (||) is a canonical transformation since it leaves invariant the 
canonical form of commutation relations. 

The vacuum expectation value of Eq.(^|) gives (|0) denotes the vacuum for the free field (p): 

<f>f(x) ee {0\i/> f (x)\0) = /+ (0| [A^d) * jW](x)] |0) . (9) 

Notice that the order parameter <f>f (x) is of full quantum nature: the classical solution is obtained by means of the 
classical or Born approximation, which consists in taking (0 1 J7" ["0 ] 1 0) = J[4>*\i i.e. neglecting all contractions of the 
physical fields. In this limit, <Pci{x) = lim?i >o & {x) is the solution of the classical Euler-Lagrange equation: 

A(d)cl> f cl (x) = J[4> f cl ](x). (10) 

Beyond the classical level, in general, the form of this equation changes. The Yang-Feldman equation ^) describes 
not only the equations for the order parameter, i.e. Eq.([lO|), but also, at higher orders in h, the dynamics of one or 
more quantum physical particles in the potential generated by the macroscopic object $ (x) [^3|. Typical examples 
of interest include for instance a scattering of quasi-electrons on the Abrikosov vortices in type-II superconductors 
or scattering of second sound waves (thermal phonons) on the vortices in superfluid 4 He. 

In refs. JT3|]l4fl , it is shown that the class of solutions of Eq.(||) which leads to non-trivial (i.e. carrying a non-zero 
topological charge) solutions of Eq. (^) , are those which have some sort of singularity with respect to Fourier transform. 
These can be either divergent singularities or topological singularities. The first are associated to a divergence of f{x) 
for |x| = oo, at least in some direction. The second means that f(x) is not single- valued, i.e. it is path dependent. 
In both cases, the macroscopic object described by the order parameter, will carry a non-zero topological charge. 

It is also interesting to consider the boson transformation at level of states [p3| . For this purpose let us write the 
generator of the field shift (|J), as 

<p f (x) = e~ w ip(x) e w = <p(x) + f{x) 

V = - [ f(y)d^(y)da». (11) 



with da^ = n^da where a is a flat space-like surface and n M its normal, both x and y belong to a. The action of V 
on the vacuum defines a coherent state |/): 

(fMx)\f) = f(x), \f)=e^\0). (12) 

The classical soliton solution is then obtained by taking the coherent -state expectation value of the Heisenberg field 
ip in the K — > limit, i.e. 

hm(/|^(x)|/> = Jim (0|^(x)|0> = <f> f cl (x) . (13) 

Although I? in (|ll]) is a precise analog of the Weyl operator of the quantum mechanics defining the canonical coherent 
states |g4|, the states |/) defined by ( fL2"|) are not the usual coherent states. Indeed, due to the fact that f(x) is 
not Fourier transformable, |/) are not eigenstates of the annihilation operator a(k) albeit they still do saturate the 
Heisenberg uncertainty relations. A discussion of the coherent states corresponding to the soliton solutions can be 
found, for example, in refs. |25|j2^]. Connection of QFT coherent states with Haag's theorem is discussed in p7|. 



B. The Haag expansion in the Closed— Time Path formalism 

Whilst the Yang-Feldman equations are quite involved and do not usually allow to proceed beyond few iterations, the 
alternative Heisenberg equations seems to be more versatile. This may be seen by considering the bare Heisenberg 
equations of the motion: 

^ B {x)=i[H,^ B {x)] (14) 
tl B {x) =i[H,n B (x)}, (15) 
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where lis is the bare momentum conjugate to the bare field ips and H is the full Hamiltonian in the Heisenberg picture 
(H = J d 3 xH(x)). The formal solution of (|l^ )-(|l5|) is well known. Assuming that the Heisenberg and interaction 
pictures coincide at some time U, we have 



xPb(x) =Z 1 J 2 K-\t) Vi„(x) A(t) 
U B (x) = Z^ 2 A-\t) Ikn(x)A(t) 
A(t) = e^- 4 ^-. e -*C*-**)« 

A(^) A~ 1 (ti) = U(t 2 ;ti) 



T 



ti 



(16) 
(17) 
(18) 

(19) 



T is the usual time ordering and Z^,, Zjj are the wave-function renormalizations (usually II oc ?/>, and so = Zn)- 
Due to the fact that both ips and V'm satisfy the canonical equal-time commutation relations, the solution ips in 
( |l6| ) must be understood in a weak sense, i.e. valid for each matrix element separately. If not, we would obtain 
the canonical commutator between ipB and TIb being equal to iZ^,^ 3 (x — y) and thus canonical quantization would 
require that Z$ = 1. On the other hand, non-perturbative considerations (e.g. the Kallen-Lehmann representation) 
require Z^ < 1. The solution of this problem is well known the Hilbert 

spaces for ips and ipi n are different 

(unitarily non-equivalent), and the wave function renormalizations Z<$ and/or Zu are then "indicators" of how much 
the unitarity is violated. This conclusion is usually referred as Haag's theorem [^1 31 . 



Let us also mention that the interaction picture evolution operator U(t 2 ,ti) alternatively reads 132 



U(t2,t!)=T* 



exp(i / d 4 x C{ n {x)) 



(20) 



where the symbol T* is called T* (or covariantjproduct^ [ p^b]j3^j33]| and C\ n is the interacting part of the density of 
the Lagrangian in the interaction picture. Eq.(E0[) is valid even when the derivatives of fields are present in C 1 (and 
thus H 1 ^ -C 1 ). 



Eq.(|16|) can be recast into a more useful form, namely 



^B{x)=Zy 3 U{U;t)i> in {x)U-\U\t) = zl /2 U(U;t f )U(t f ;t)^ m (x)U(t;t t ) 



Zl /2 T C 



tpi n (x) exp I — i / d 4 xTii n (x) 



= 4 /2t c 



ipin(x) exp ( i / d i x£ I in (x) 



(21) 



Here C denotes a closed-time (Schwinger) contour, running from ti to a later time tf, t < tf and back again (see 
Fig.l). Similarly, Tq denotes the corresponding time-path ordering symbol (analogously for the Tq ordering). In the 
limit ti — > — oo, we get that ipi n turns out to be the usual in-(or asymptotic) field. As the time tf is by construction 
arbitrary, it is useful, from a technical point of view, to set tf = +oo. Eq.(|2l|) may be viewed as the Haag expansion 
of the Heisenberg field -0s- 



Generalization of Eq.(21) to the case of more fields is straightforward. So, for instance, for the time ordered product 



of n Heisenberg fields we may write 



T[lp B (xi) ■ ■ ■ 1pB(x n )] = Z" /2 Tc 



4>i n (xi) . . . tpi n {x n ) exp ( i I d 4 xC{ n (x) 



(22) 



When we consider vacuum expectation value of Eq.(|22|), we have at our disposal two equivalent representations for 
the T — Green's functions, namely 



2 The T* product is defined in such a way that for fields in the interaction picture it is simply the T product with all 
the derivatives pulled out of the T-ordering symbol. Evidently, for free fields without derivatives: T*[tpi„(xi) . . . ipin(.x n )] = 

T[lp in (xi) . . . 1pm(Xn)]. 
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FIG. 1. The Closed Time Path 



(0\T(^ B (x 1 )...^ B (xn))\0) = Z , / 2 (Q\T^ ^ n { Xl )...^ m {x 2 ) exp(ij^d 4 x£f n (x)\ 



10) 



(23) 



_ ~n/2 

~ Z 4> ' 



(0\T* 


lpin(xi) ■ 


■ *Pin{x n ) exp (i f™ x d A xC{ n {x)} 


10) 


(0|T* 




10) 



(24) 



Eq.(p4[) is the well known Gell- Mann-Low formula for Green's functions. Note that the latter is true only for the 
vacuum expectation values, due to the stability of the vacuum |0): 



exp i 



with a being an eigenvalue (basically a phase factor) of the interaction-picture evolution operator which corresponds 
to |0). In more general situations where expectation values are taken with respect to a state \ip) ^ |0), Green's 
function cannot be recast in the form (|24|) and the T<S> prescription is obligatory. So for mixed states, \ip) — > p, 
(il>\...\i/>)^Tr(p. ..) = (•••) and 



{T{P[M)) = (T*c 



P r [ipin] exp 



d 4 x 



(25) 



where P r [. . .] is an arbitrary (generally composite) polynomial in ipg, and the subscript r suggests that the correspond- 
ing renormalization factors are included. Eq.(p5|) will be of a fundamental importance in our following considerations. 

An important special case is that of a system in thcrmodynamical equilibrium. Then the statistical properties of the 
system are described by the canonical density matrix (for simplicity we omit from our consideration grand-canonical 
ensembles). As p cx e~@ H , the density matrix is basically a generator of the (imaginary) time translations. Using 
Eq.(|l8|) we may then write 



-0H 



U(ti-i/3,ti). 



From this it is evident that one may substitute the full density matrix with the density matrix for the corresponding 
free system provided one adds to the path C a vertical part running from ti — i(3 to ij (see Fig. 2). Advantage of 
this rather formal step is that the free density matrix is Gaussian and correspondingly (Tc[ipi n (xi) ■ ■ ■ tpi n {x n )]) is 
for n odd and a symmetrized product of the (free) two-point Green's function if n is even. This is nothing but the 
thermal (or thermodynamical) Wick's theorem - basis for a perturbation calculus in QFT at finite temperature. We 
shall elaborate more on this point in Section IV. 



C. Heisenberg operator in the presence of defects 



Having obtained a compact expression for the Heisenberg field in terms of the interacting fields, we can select the 
initial-time conditions corresponding to a particular physical situation. In other words (see Section II. A), we select 



G 



FIG. 2. The Thermal Path 



the set of asymptotic fields describing the particle content of the theory in question. For instance, in the case of 
a type-II superconductor we must choose either quasi-electrons, if we are above the critical temperature (normal 
phase), or bogolons if we are bellow the critical temperature (superconducting phase) p4j. It is well known that the 
new vacuum state in the superconducting phase is an (homogeneous) condensate of Cooper's pairs (BCS state). 

Similarly, for paramagnetic-ferromagnetic phase transition in the Heisenberg model we choose particles (atoms) with 
spin i above the Curie temperature (paramagnetic phase) whilst we choose magnons if we are below the Curie 
temperature (ferromagnetic phase). The vacuum in the ferromagnetic phase is then a state where all particles have 
spins aligned in the same direction - an (homogeneous) spin (or Bloch) coherent state. 

Let us also mention the important case of the QCD phase transition. Here, in the confined phase (i.e. at energy scale 
<S Aqcd) one usually chooses as asymptotic states hadrons and mesons (i.e. quark bound states). In the deconfincd 
phase, i.e. at temperatures of order ~ 0.2GeV ~ 10 12 K where a quark gluon plasma starts to form and the chiral 
symmetry is restored, one could, in principle, directly work with quarks and gluons as relevant asymptotic fields. 
In this case, the structure of the confined phase vacuum is still far from being understood. The latter is sometimes 
identified with the (homogeneous) chiral condensate or with "color" BCS state, etc. 

The choice of the initial-time data is therefore an important theoretical probe for the description of phase transitions 
and their underlying microscopic mechanism. If a phase transition is sufficiently rapid it may happen that stable 
topological defects are created; examples include vortices in rotating superfluid 4 He below the A point, a rich variety 
of defects in 3 He, quantum magnetic flux tubes (Abrikosov vortices) in type-II superconductors, disclination lines 
and other defects in liquid nematic crystals, droplets of an unbroken phase in quark gluon plasma, etc. As we have 
already mentioned in Section II. A, such topological defects can be generated with a special choice of the initial-time 
data by means of the boson transformation with a Fourier non-transformable shift function /. Loosely speaking, 
the emergence of topological defects may be seen in the corresponding vacuum state for the asymptotic fields which 
should be a "suitably chosen" inhomogeneous coherent state. 

We can summarize our strategy as follows (see Fig. 3): For a given dynamics, the CTP formulation gives a closed 
(functional) expression for the Heisenberg field operator(s) -0 m terms of the asymptotic (physical) fields. We then 
use the boson transformation to introduce the shift function / controlling the choice of the Hilbert space. 

The next step is to consider the order parameter, i.e. <fif = (0|?/>^|0). By taking the classical limit (j> cl = limj i _ + o <f>f , we 
determine the form of the shift function, say /, corresponding to a particular (classical) soliton solution. This shift 
function is then used to obtain the Heisenberg field operator in the chosen soliton sector: ip* . At this point there are 
various possibilities: 

- calculate quantum corrections to the order parameter, by taking higher orders in the h expansion of (Olf/^lO); 

- study finite temperature effects on the order parameter, by considering (t/jf)^, where (...) p stands for thermal average; 

- calculate Green's functions in the presence of defects, both at zero and finite temperature, as (T[i/v i/v ]), etc. 

A more ambitious task is the study of the non-equilibrium properties of QFT systems containing defects, e.g. consider 
phase transitions. This will be the object of our future work. 
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Heisenberg field Boson transf. v.e.v. 

Win] * ^ = </#» + /] * ^=(01^10) 



classical limit 



choice of f 

Heisenberg field Jjf f J c l assica l so Uton 

in presence of defects r J Yd classical suiuon 



T = ► (0|^|0) , (0\i)flpf\0) , ... 

7V0 ► (ipf) , (^V~)/3 , - 

p + e-^\ ► (tf) , (^V) , ... 

FIG. 3. Iterative procedure for computations 



III. KINKS IN TWO-DIMENSIONAL \ip 4 THEORY AT T = 

In this Section we apply the formal considerations developed above to a specific model, namely the two-dimensional 
A-0 4 theory at zero temperature. It is well known that Goldstone's sector of the Aip 4 theory in (1 + 1) dimensions 
possesses at a classical level a certain class of analytical extended solutions, namely the kink (antikink) solutions jiB). 
In the following, we derive the classical kink solutions at T = and construct the Heisenberg operators in presence 
of kinks. As a further non-trivial step, we extend in the next Section this analysis to finite temperature. 

A. The Heisenberg operator in the presence of a kink 

Let us consider the case of a 1 + 1 dimensional scalar system of the hermitian field tp with quartic interaction term. 
We have the following bare Lagrangian (throughout we adopt the Minkowski metric with signature (H — )): 

= \{d^) 2 -\^ 2 -- A ^ + C ctl (26) 
where C c t is the counterterm Lagrangian: 

Cot = \(Zip - 1) {d^f - \ W -\{Zx- 1)^ 4 • 

Here we have introduced in the usual way a renormalized field, mass and coupling as: ips = Z^/ 2 ip, Z^,^ 2 B = /i 2 + S/.1 2 
and Xb = Z^ 2 Z\X. The corresponding renormalized Heisenberg equation of the motion for the field operator ip is 
(see e.g. g§) 

(d 2 +^(x) = -\{^(x)}. (27) 

The squared bracket on the RHS of (|27]) denotes a renormalized composite operator (or Zimmerman normal ordering) 
corresponding to an unrenormalized operator ip 3 . 

In the case when both /i 2 < and (Aq < (i.e. in the Goldstone phase) it is well known that ( pif ) admits at a classical 
level kink solutions pq , ^ , 06| . We shall omit from the next discussion the Wigner phase (i.e. when > 0) since it 
does not enjoy the kink solutions. Let us define the following renormalized quantities 
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ip(x) = v + p{x) 



with v = 



-— , -2p = m 



g = V2X. 



(28) 



The parameter v represents a tree-level expectation value of ip (i.e. the true classical minimum of the potential in 
(pq)). As a consequence of the parameterization (pq) we obtain 



l/„s9 loo Q A 1 o m 



S 



8ff 2 



(29) 



where 



£ct = ~ - l)(d, P f \ \{Z X 1)| to 2 - 5m 2 J p 2 - ~ ,g(Z A - l)mp 3 ~\{Z X - l)p 4 
- i ™ (<5to 2 - (Z A - 1)to 2 ) p + i ^ f <5m 2 - 1 (Z A - l)m 5 

Here we have used that 8m 2 = —28pL 2 . The form ( p9| ) of the Lagrangian is particularly important in doing perturbative 
theory. We may recall that both m and g in C are clearly (by construction) arbitrary parameters and no physical 
observable should depend on their choice. But, in practice, when one deals with perturbative calculations, choice 
of one parameterization over another can save some labor and possibly cultivate on the intermediate stage a better 
physical intuition. Here, and in the work to follow, we explicitly choose the zero-momentum mass renormalization 
prescription: i.e. Ei?[p 2 = 0] = —to; with E/{ being the renormalized self-energy. Similarly, we chose the zero- 



momenta coupling renormalization condition: T 



(4) 



(4) 



[0,0,0,0] = -A; with T R 



being the renormalized 4-point proper 



vertex function. Both m and A are chosen to be the tree level mass and coupling, respectively. The virtue of this 
prescription is that we keep from the very beginning the interpretation of v as the tree-level vacuum expectation 
value of ip. In addition, the zero-momentum renormalization prescription can be easily recast into conditions on the 
renormalized effective potential which is the most natural tool for a treatment of phase transitions |50|,[37j (this is 
particularly useful when comparing with existing results which are predominantly computed in the effective potential 
framework). The physical, but mathematically more involved on-shell renormalization may be obtained via a finite 
renormalization procedure ]30| ] . Let us also remark that in 2D case the A ^ theory is super-renormalizable and so 
both Z\ and Z$ are finite to each order in the perturbative expansion, while the divergence present in <5to 2 is solely 
due to the tadpole diagram which is here only logarithmically divergent. 



The renormalized Heisenberg equation for the field p is: 



(d 2 + m 2 )p(x) = ~\mg[p 2 (x)\ - l -g 2 [p\x) 



(30) 



Note that if p 2 < then m 2 > 0. The asymptotic (t — * — oo) field can be then identified with the free massive field 
satisfying 

(d 2 +m 2 )p ln (x) =0, 



(31) 



with to being the aforementioned tree level renormalized mass. Setting in Eq.(|l8|) Hf n = \p 2 in + \p' % n — \ m2 p 2 n then 
the interaction Lagrangian C 1 (entering Eq.pO)) is 



£ J [p] = CM - ~ g 2 p 4 - \ mgp 3 + ^ , 



and the dynamical map for the field ip has the form 



^{x) =v + T* c 



p in (x)exp\i / (fy^lpin] 



c 



(32) 



Note that the constant factor in C 1 automatically cancels during the contour integratiorj^]. Let us now consider the 
boson transformation: 



This fact, among others, shows that in the CTP formalism the UV divergence in the energy density of the vacuum is 
automatically cancelled. 
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Pin(x) -> Pin{x) + f(x) 

(d 2 + m 2 )f(x) = 0. 
As a result we get the new Haag expansion for the field (x) 



(33) 
(34) 



ip f ( x ) = v + Tl 



(/.,„[./■) f f(x))exp{i 1^ ■!-!,£' ,,,„[</> ■ f(z)] 



iSJ(x) 



/(*) 



exp<j i I d 2 yC I 



iSJ(y) 



Tc expi / d yJ{y)p m (y) 



(35) 



,7=0 



where we have introduced the c-number source J in order to perform some formal manipulations. By use of the 
(operatorial) Wick's theorem pO] , we get 



ip*(x) = v + 



iSJ(x) 



expti I d 2 yC I 



c [iSJ{y) 



■M 



x :expi / d 2 y J(y)p ln (y) : exp 



~ f d 2 yd 2 zJ(y)A c (y;z)J(z) 



(36) 



where Ac{x;y) — (0 1 Tgipm (%) Pin (ll) ) 1 0) and |0) is the vacuum for the pi„ field. Once the function / is "properly" 
chosen (see below), Eq.(|36[) represents a convenient representation for the Heisenberg operator in presence of defects 
which can be used for further analysis |la,M . 



In order to determine the function / leading to kink solutions, let us now consider the vacuum expectation value of 
the Heisenberg field ip' . Here the normal ordered term drops, and we get ((. . .) = (0| . . . |0)) 



(ip f (x)) = v + 



iSJ(x) 



expii I d 2 y£ I 



c 



*SJ(y) 



f(y) 



x exp 



d 2 yd 2 z J(y)A c (y;z)J(z) 



.7=0 



By use of the relation 



iSJ 



G[J] = G 



iSK 



F[K]e^ KJ 



(37) 



(38) 



K=0 



we obtain 



tyi* (x)) — v + exp 



d yd zA c (y;z)- 



i8K(y) iSK(z) 

x (K(x) + f(x)) exp \ i I d 2 y [C T [K(v) + f(y)] + K(y)J(y)] 



c 



(39) 



K=J=0 



We now perform a change of variables K{x) — > K{x) + f(x) and set to zero the term with J (there arc no derivatives 
with respect to it). As a result we obtain 



{tjjf (x)) = v + exp 



-- j^d 2 yd 2 zA c {y;z)- 



8K{y) iSK(z) 



K{x)B[K] 



K=f 



with 



B[K]=expii / d'yC^Kiy)] 



c 



(40) 



(41) 
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We can thus express Eq. ( ptcf ) as a sum of three terms [j] 

(^(x))=v + C[K](x)\ K=f + D[K](x)\ K=f 



C[K](x) 



d 2 yA c (x;y) 



c 



5K(y) 



exp 



d zd yA c (z;y) 



5K{z)5K{y) 



B[K], 



D[K](x) = K(x) exp [\ £ d 2 xd 2 yA c (x; y) } . 



(42) 
(43) 
(44) 



Let us observe that C[K — 0](x) is independent on x since the inhomogeneity in the order parameter is controlled by 
the x dependence in K(x) (or precisely in f(x)). As a result we see that C[0] describes quantum corrections to the 
tree level homogeneous density of condensation v (note D[K = 0] = 0). Introducing the notation: 

C[K](x) = C[K](x) - C[0] , v = v + C[0], (45) 

with w being the renormalized order parameter, we arrive at the following alternative form for the order parameter: 

(^(x)) = v + C[K]{x)\ K=f + D[K](x)\ K=f . (46) 

Note that both D[K] and C[K) now vanish when we set K — 0, which is equivalent to considering the homogeneous 
condensation only. In the literature, C[0] is often denoted as Sv P8| . 



B. The kink solution in the Born approximation 



So far, all result obtained were of a full quantum nature. We now deal with the Born, or classical, approximation of 
Eq. ( (46[) and for this purpose we reintroduce h. Each propagator is then proportional to h whilst B has in the exponent 
the factor ihr 1 . The Born approximation means than only terms of order h° in ( f46|) must be taken into account. In 
order to better understand the classical limit, let us view the held theory as a classical physicist does: one would then 
talk about a frequency loq and a wave vector k of the field rather than about mass m and momentum p. In QFT, the 
particle mass and momentum are obtained by multiplying ljq and k by h (we do set c = 1). The units of A = g 2 /2 
are [i? 1 ^ -3 ] where E means units of energy and I means units of length, so [A] = [ftr 1 ^ 2 ]. Note that the latter also 
implies a super-renormalizability of the theory. Similarly, the Fourier transform is expressed in terms of the product 
kx and not px. These comments, even if trivial, will be of a particular importance in the finite temperature case 
where only a carefully performed classical limit can provide consistent results. For instance, in massless theories, a 
correct classical limit reveals the breakdown of classical perturbative expansion (infrared catastrophe) and then some 
resummation is required - the so called hard thermal loops (HTL) resummation |38|,^9| . 

In rewriting some of the previous expressions making the h dependence explicit, we recall that the counterterm 
Lagrangian C c t is of order 0(h) and so in the expansion of C[f] and D[f] terms one always gets loop contributions. 
It is obvious that in doing classical limit such contributions must decouple and so in the following reasonings we may 
neglect £ ct . We have n 



^ d 2 y A c (x; v) 5K ^ ex P \ ha \ ex P 
D[K](x, h) = K(x) exp [ha] exp 



(47) 
(48) 



'We use the identity 



exp 



- ^2 Aijd Xi d Xj x k B{xi)= ( x k + ^2 A k jd Xj J B(xi + ^ Ai 3 d Xj )l, 

ij J \ j I i 



where Xk —> K(x) and B(x) — » 
5 After factorizing out h, we shall in the following consider Ac(. . . ) oc h°. 
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with 



a =\J d 2 zd 2 yA c (z;y) 



S 2 



b= -r / d 2 i 

c 



5K{z)8K{y) 
[ ^K\z) + lo o9 K\z) 



(49) 
(50) 



Keeping only the finite terms in the fi — > limit we get[] 

C[K](x, h -> 0) = y d 2 yA c {x; y) j^r^ Resr i= o fexp [fia] exp ^6 



= t d 2 yA c (x;y) — V ,, 1 | a 

/oo (5 00 1 

d 2 y Gr(x, y) — — V a"6 



(51) 



where iGr(x, y) — 0(xq — yo) A(x; y) is the (position space) retarded Green's function of the free theory and A(x; y) 
is the Pauli- Jordan function: 



A(x o ,xi;0,0) = (0|[p 

in (^) 5 Pin (0)]|0) 

The D term gives 



|A^(fc 2 -c 2 ) £ (fc )e-^ = -^ 



8(x - |xi|) J (w ^x 2 , - x\) . (52) 



-D(x, h —>■ 0) = if (x) Res^o ( — exp [Tia] exp 



1 



OO 



n=0 



In Appendix A we show that -D[if, H — ► 0] = K. The final result is then 

S 



/oo A 1 



E p »^K a 



where 



Pi(x) = K{x) 

Pn(x) 



1 



d 2 yG R (x,y) a r ' 



^u g K 2 (y) + \ 9 2 K 3 (y) 



n > 2. 



[{n-2W 

In Appendix C we show (using mathematical induction) that the following recurrence relation holds: 



Pn(x) 



d 2 yG R (x,y) 



n > 2, 



i-^j—n i-\-j-\-k=n+l 

where i, j, k = 1, 2 . . . . The recurrence condition ( |56|) may be "diagonalized" in the following way 



(53) 



(54) 



(55) 



(56) 



6 Throughout this paper h — > is a short-hand notation for ft — > O4 
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h(x) 
h(x) 



K(x) 

/CO o 
d 2 y-G R (x,y) E QiMQM 
-OO ZV A 1 A—I 



i+j=2 



d 2 yG R (x,y) 



E QMQM 

i+j=n 



\g 2 E QMQMQkiv) 



i-\-j-\-k—n 



■ n > 3. 



with 



Qn — P%Pn+l + Pn 



(57) 



(58) 



Note that both ^» anc ^ Si^i l ea, d to the same result^and so it does not really matter which recurrence 

equation will be solved. The formulation of the problem in terms of recurrence conditions ( |57| ) has however one 
crucial advantage, namely Eq.(|57j) belongs to the class of the so called functional equations of Cauchy-Marley's type 
pO| , p| where it is known that the fundamental solution cannot be expressed (apart from a very narrow class of 
kernels) in terms of elementary functions. This means that we cannot resolve ( |57| ) in terms of general K (or /), i.e. 
we cannot find the Backlund transformation between the solutions of Eq. ( |3~I| ) and those of Eq. (|64|) - the classical Euler- 
Lagrange equation of motion (see below). Nevertheless we can obtain all the analytical solutions of Eq.(|(34|), namely 
the analytical kinks: This is done once we realize that the convolution of the 2D retarded Green's function Gr(x] 
with an exponential is proportional to the very same exponential (actually the exponential is the only function having 
this property - see Appendix C). So we choose Q n {x) cx (Qi(x)) n with Qi(x) = fo(x) being an exponential, Fourier 
non-transformable solution of the Klein-Gordon equation (B4j) . We obtain: 



Qn(x) = A n f£(x) = A n e ± ^ n ^~ x ^ , 



(59) 



where 7 = (1 



(u will be later interpreted as a velocity of a kink). Plugging this form into the recurrence 



relation (57) and using the result from Appendix C we arrive at the following equation for the factor A n 



(ri 



l_ J+j=n i+j+k=n 



AiAjAk 



(60) 



This is a trivial version of Cauchy-Marley's equation which has the only (non-zero) fundamental solution, namely 
A n cx (Ai) n . Using the standard identities: 



i-\-j=n 



E l = \(n-l)(n-2), 

i-\-j-\-k—n 

it might be easily checked that a solution of the recurrence relation reads 

s 



At. 



2v 



2v 



with c being a real constant. Thus, finally, we have 

'sfo(x) 



2v 



= vcth 



2 V 2v 



(61) 



(62) 



(63) 



Here Ln(z) = ln|z| + i argz is the usual principal value of the logarithm of z. Note that K = f = s/q. Thus, provided 
f(x) is an exponential solution of the linear equation (0), the solution ([33|) fulfils the (classical) Euler-Lagrange 
equation of motion: 



7 This might be immediately seen from the fact that both 2~2i=i P» an< ^ X^i^i Q* solve the same integral equation, namely 

1 



X(x) - f{x) 



d yG R ix,y) 



*-u 09 X\y) + ^g 2 X\y, 
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-w2/2. 



(64) 



The later is nothing but the expectation value of equation (|27| ) in the Born approximation. 

-2ve UJ °' ya , we easily obtain the standard kink solution 



For instance, if we choose fo(x) = e "otC 351 with s 

H-0 



(V'oW) =wth y7((xi - a) - x u) 



(65) 



describing a constantly moving kink of a permanent profile with a center localized at a + ux$. Note that the former 
function f(x) is the solution of the homogeneous Klein-Gordon equation (34), is not Fourier transformable and fulfils 
the initial value condition: f(xo — ► — sgn(u)oo, x\) — > 0. 



As the Lagrangian (pq) is Z2 invariant we could equally choose (0|^(a;)|0) = —v. In this case we would get 



(tpoix)} = -vth -^7((xi - a) - x Q u) 



(66) 



which is the antikink solution. Note that the antikink can be also alternatively obtained choosing fo(x) 
e+^nixi-xou) and s = _ 2t;e -"o7a j provided we keep (0\ip(x)\0) = v. 



It is quite instructive to realize what happens if we choose s positive, i.e. s — 2ve u °' ya . In this case we obtain 



(x)) = ±v cth ^7((xi - a) - x u) 



(67) 



(± sign depends on the choice of fo(x)) which is also solution of Eq.(Q). However, this solution is singular in 
x\ — a = xqu and so it does not classify as a soliton even if, for example, the corresponding energy density 



H cl {x) = \{{^{x))) 2 + l^x)) 1 ) 2 - \w 2 M(x)) 2 + \{^{x)f + \u:y 



2 2 i4 

-v LOr, cscn 



j((xi - a) - x u) 



4(1 — 

is localized near x\ — x u + a, and the corresponding total energy 

,3 



/ U! 
dx x H cl (0, Xl ) = ^7 = M ca ., 



(68) 



(69) 



is finite (M c i denotes the classical (static) kinkmass). 

We thus see that, in the classical limit, the present method gives us a constructive way for finding classes of analytical 
solutions of non-linear equations. The fact that we could not find the fundamental solution of Eq.(|34|) is intimately 
connected with the theory of Cauchy-Mar ley's functional equations. On the other hand the non-existence of the 
fundamental solution is fully consistent with the fact that (1 + 1) dimensional Xip 4 theory is not exactly solvable. 

The above solutions for the shift function / corresponding to kink solutions can be now used in Eq.(|36|) to get the 
Heisenberg operator in a given kink sector. One can thus calculate Green's functions and higher quantum corrections 
to <f>f . We show this elsewhere [ttq,f9[ . 



C. The kink solution in the small coupling limit 

We can obtain the classical limit of the order parameter (Q) using a different standpoint, namely the semiclassical 
(or WKB) approximation . This asserts that a system behaves classically if its de Broglie wavelength A^s is much 
smaller than the characteristic length scale ro of the theory: at this stage the wave properties of a system become 
unimportant. In our case the kink (antikink) de Broglie wavelength is XdB = 2irh/\p\ = 2irh/~fM c i\u\. On the other 
hand, ro in the theory can be identified with the kink (antikink) width (= I/7W0). So 

Ads < r Q =*> hX <C WqM & % 2 < w oM- ( 70 ) 
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Note that the above condition can be fulfilled in all Lorentz frames of reference apart from the ones where the kink 
(antikink) is close to or exactly at rest. This has a simple physical meaning: the semiclassical approximation is 
based on the assumption that typical momentum and distance scales are much larger than the error scales saturating 
Heisenberg's uncertainty relations. This assumption naturally fails when the typical momentum scale is close to zero 
and correspondingly for small |u| the semiclassical approximation cannot be used. 

Eq.(f70|) is fulfilled either if h — > (see previous subsection) or if A — > 0, i.e. in the heavy mass approximation. 
We study the latter case in Appendix D: there we show that the small coupling limit, leads to the analytical kink 
(antikink) solutions of the classical Euler-Lagrange equation (|64|). We stress once again that the classical soliton 
solutions obtained both in the previous and present subsections were determined as a classical limit of the full 
quantum mechanical expression for the order parameter (|46|). 

IV. KINKS IN TWO-DIMENSIONAL A i> 4 THEORY 

So far we have dealt with QFT at zero temperature. However, in many physical situations, the zero temperature 
approximation is not appropriate. This happens for example in cosmology [|46|-[49|| , astrophysics [^p|^2| and in the 
study of the quark-gluon plasma formation J5^-|55| as well as in condensed matter. More in general, in all cases 
when one studies systems near the critical temperature (and so neither low or high temperature expansions properly 
describe the basic features) the finite-temperature treatment must be carefully taken into account. 

It is important therefore to understand what happens with the topological defects when the QFT system is immersed 
in a heat bath (thermal reservoir) which is saturated at some temperature T. In this Section we aim to use our toy- 
model system to address this question and demonstrate the basic features of our approach in the finite-temperature 
case. Most of the results here obtained retain their validity when more realistic (higher dimensional) systems are 
considered. 

As already discussed in Section II. B, in thermal equilibrium, the most convenient choice for the time path is the one in 
Fig. 2, to which we refer as the thermal path. The crucial observation at finite temperature is that the operatorial Wick's 
theorem still holds (see e.g. |5(|) and consequently Eq.(|3^) retains its validity provided the following substitutions are 
performed: 

A c (x;y) = (0\T c (p in (x)p in (y))\0) — > A c (x;y,T) — (T c (p in (x)p in (y))) p 

— N(...), 

where (. . .) p = Tr (e~P H . . . ) /Tr (e _/3ff ) and (3 — 1/T. The thermal normal ordering N(. . . ) is defined is such a way 
f36f that (N(. . . ))p — 0, the dots stands for a product of T = free fields. This is of a great importance as all the 
formal considerations developed in Section III. A go through also for finite T. 

A. The kink solution at finite T in the classical approximation 

The classical approximation has proved to be a useful tool for study of low-energy and/or high-temperature properties 
of quantum fields at finite temperature. Examples include non-perturbative computations of the Chern-Simons 
diffusion rate JstJ , sphaleron decay [Bp or plasmon properties of hot non-abelian gauge theories |59| ; there is also a 
perturbative treatment via the classical Feynman diagrams |]60|| which however will not be followed here. 

At nonzero temperature the question of h appearance is more delicate than in the zero-temperature case. The whole 
complication is hidden in the thermal propagator Ac(x;y,T). Whilst at T = the latter is directly proportional to 
h, at finite T the situation is very different. Although this generic observation appears simple, there has been in the 
past confusion in the literature about this point which was understood properly only recently p§| , |6l| . To understand 
the complications involved let us make H explicit. The free thermal propagator in spectral or Mills's representation 
||,|9| then reads 

A c (x;y,T) = h J ^e- ik ^p(k)[e c (x -y ) + f b (hk /T)} 

= A c {x;y) + A T c {x;y) : (71) 
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with 



A c (x; y) = hj 0- e - ik ^p(k) [0c(x o - yo) - 0(-k )] 



(72) 



where the spectral density p(k) = (2n)s(ko) 8{k 2 — ujq) with e(fco) = 6(k ) — 8(—k ). The contour step function 
8c(xo — Uo) is 1 if yo precedes xq along the contour C. The Bose-Einstein distribution fb[x) — (e x — l)~ . It might 
be directly checked that we obtain the usual elements of the (free) thermal propagator in the so called Keldysh- 
Schwingcr formalism p3,B9|, i.e.: 



A n (x;y,T) 



d 2 k 



-ik(x-y) 



2ir6(k 2 -u; 2 )f b (h\ko\/T) 



(2tt)2 

A 21 (x;y,T) = h J ^ e ~ ik ^ {9(k ) + f b (h\k Q \/T)} S(k 2 - oj 2 

A 22 (x;y,T) = (A 11 (x;y,T))* 
A 12 (x;y,T) = (A 21 (x;y,T))*, 



(73) 



where "*" denotes complex conjugation and the particle mass m = Hloq. We remark that the thermal part of Aq(x; T) 
is identical for all matrix elements and that in ([73]) it appears fb(h\ko\/T) and not fb(Hko/T). We should also emphasize 
that k in the integration is a wave vector - a reciprocal length - and not a momentum. 

Now, due to the mentioned analogy with the T = situation, we may immediately write for the order parameter 

{^{x))p = v + C[K](x-T)\ K=f + D[K](x;T)\ K=f , (74) 



where we took the thermal average of the expression analogous to the one in Eq.(|36|), but with the normal ordering and 
the propagator replaced with their thermal counterparts. Both C[K\{x;T) and D[K](x;T) entering Eq.(|74|) coincide 
with their zero temperature counterparts provided one uses the thermal propagator instead of Ac(x,y). 

Let us consider the classical limit of Eq. ([74]) . The D[f](x;T) term then gives: 

1 (h 



D[f](x;T,h^0) = f(x) Res h=0 
x exp 

= f(x)Res h=Q ^2 



2h 



d 2 z 



n.rn— 



9 1k< 

4 



K=f 



K*(z)+Lu gK 3 {z) 

X - f d 2 zd 2 yA c (z;y,T) 
J c 



8K{z)5K{y) 



(75) 



K=f 



Note that if Ap(i; y) = we recover the result (|53| ) of the previous section. Using the result of Appendix A we may 
directly write 



D[f] (x; T, h^0) = f(x) Res h=a ]T 



SnO 8 



mO 



fix). 



(76) 



For the C[f](x; T) term we get, in the classical limit: 
C[f](x;T,h^0) 



= J^yAdx^—Res^Y, 



— I d 2 zd 2 w Ac(z;w,T) 
2 Jc 



SK(z)5K(w) 



(77) 



K=f 



There is no simple way how to evaluate ( f77j ) without performing an explicit Laurent's expansion of Ac{x 7 T) around 
h = 0. Using the Bernoulli expansion: 



x = y R El 
eP - 1 ^ a a\ 



\x\ < 2tt. 



a=0 
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(B a are Bernoulli's numbers) we may Laurent expand fb as 

, , 1 T 1 hk 1 fhk Q \ 2 , . 

The key observation is that the foregoing series converges only for h\ko\ < 2-kT. The leading term in ( £q ) gives the 
classical thermal part of the propagatoijj, and the classical approximation is then equivalent to taking the leading 
term in the Laurent expansion ([78]). The higher quantum corrections are due to the higher terms in the expansion, 
but for large |fco| the expansion does not work, i.e. an expansion in h is unwarranted. Of course, for 7i|fco| 3> T the 
distribution fb is exponentially small (Wien's distribution law) and it is dropped in comparison to the (zero point) 
first term in the integral (J7l|), which returns the usual T — approach from the previous Section. 



In order to get some quantitative results let us concentrate on the high temperature case (low temperature case is more 
involved and will be treated elsewhere). By high temperature we mean the temperature at which the Rayleigh-Jeans's 
sector of the thermal propagator (|7l]) approximates by itself sufficiently well the thermal part of the propagator. This 
should be taken with a grain of salt as the high temperature we intend to use should be still bellow the temperature 
threshold above which the heat-bath quasiparticles become instable due to the thermal fluctuations. 

Because we are in the high temperature regime we take cut off in the ko integration: \ko\ -C Ai » T/H. Due to 6- 
function in the thermal propagator we must also impose cut off on k\ integration: |fci| <C A2 ~ T/h. We tacitly assume 
that the domain of integration from T/h to 00 will contribute modestly due to the exponentially small contribution 
to the thermal propagator. 

Using (|5l|), ( |63| ) and ( |77| ) we can now evaluate C[f](x; T,H—+ 0). Indeed 

5 



C[f] (x; T, h 0) = exp [c] li d l y Gr(x, y) 



SK(y) 



r 4v 2 

exp [c] <^ -K(x) -2v- —— — 

l_ K (x) — 2v 



V - a n b n+1 1 

^ n!(n + l)! f 

n=0 V ; ) 



K=f 



(79) 

K=f 



where 



Note that there is no h in c in the large T approximation. In other words, the high temperature properties are injected 
into the order parameter calculation in a purely classical way. This might be compared with the usual observation 
that at large T thermal fluctuations dominate over quantum ones. Thus in our case we get (we omit h in as it 
cancels anyway in the large T limit 



(81) 

K=f 



Applying the expansion 



1 



= A' 1 - A~ 1 BA~ 1 + A^BA^BA- 1 



A + B 

'SK(y) 



and identifying A = K{x) — 2v and B = J c d 2 y A^(x; y) S J, . we get 



8 We can call it Rayleigh-Jeans's sector of the thermal propagator as the corresponding distribution function f{uj) = T/lu is 
nothing but Rayleigh-Jeans's distribution law. 
9 To derive Eq.@ we use the relation: exp fe^ T,'>, i). ) F(x k ) = F{x k +2£,. A k jd Xj )l. 
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riru t i (1 , n, •> 1 , A^(0) [A£(0)] 2 
C[/J(x; T, n -»• 0) = -f{x) ~2v-4v — — — — + — — — - + 3 



K{x) - 2v (K(x) - 2vf (K(x) - 2v) 5 



~ /l j /(a;)-2t; £j (K (x) - 2v)^ 1 j ' 



K=f 

(82) 



The series = 1 + S^Li z n (2n — 1)!! = X)^o cr ™ z ™ ^ s clearly divergent for all non-trivial values of 2 = 

A^ 1 (0)/(.K"(a;) — 2v) 2 . Thus can be at best understood as an asymptotic series. Question arises then which 

function has S(z) as its asymptotic expansion. To determine this, let us assume initially that zeC and perform the 
Borel transform Bs(z) of by |3] 



B,( 2 ) = S B„ 2 » sE ^= (83, 

n— n— v 



Consequently the Borel sum of E(z) is 



^e-^( te )=-^e-V- Er fc(- 7 i=) ) (84) 



where Erfc(. . . ) is the complementary error function (see e.g., |65J]). Relation ( p4j ) is true for any z € C with arg(z) ^= 0, 
and so in such a region S(z) is Borel summable with the unique asymptotic function 



e' 1 ^- Erfc (y^J ~ , arg(z) ^ . (85) 



In the marginal case when arg(z) = 0_ or arg(z) = 0+ we have 



e 1/2z Erfc ( i-^= ) - S(z) or i\l—e~ 1/z - U/^- e _1/2z Erfc ( ) ~ S(z) , (86) 



2z \ y/2z J V z \ 2z \ y/2z_ 

respectively, depending which value of \J — 1 one accepts along the branch cut. From the uniqueness of the asymptotic 



expansion (p5j) follows that when z is real, then the functions in (86) are the only ones which have £(2) as their 
asymptotic expansion. This dichotomy is inherent in theory of asymptotic expansions fl64| and fortunately is not 
counterproductive here. Both functions differ by an exponent which rapidly vanishes when z — > and accordingly 
they belong to the same equivalence class of asymptotic functions (see e.g., Poincare criterion [^4[ ). A brief inspection 
of (pq) reveals an important feature of the class functions, namely that they are complex even when z is real. 

Applying now the fact that in (1 + 1) dimensions the thermal tadpole A^ L (0) = T/2ojq = T/2gv, we see that there is 
no high T solution for the thermal kink as the order parameter turns out to be a complex number. This result is fully 
compatible with the absence of phase transitions in (1 + 1) dimensions p^j67| . In fact, the emergence of the complex 
valued order parameter can be attributed (similarly as in the case when the effective action machinery is in use) to 
spontaneous fluctuations in the order parameter. The latter behavior is well known, for instance, from Ising-type 
models f(38|| or from numerical simulation in (1 + 1) dimensional 4 theories [ ^0[ . 

It is interesting to compare the above result with the one contained in Ref. there a tree approximation at finite 
temperature is used to get a "classical" kink solution which is of the same form as the zero-temperature one but 
with thermal parameters m(T) and A(T) in the place of zero-temperature ones. From our discussion it is clear that 
the tree approximation is not useful in extracting classical results at finite temperature as the loopwise expansion 
is not anymore an ex pan sion in h. In fact, a resummation of infinitely many (thermal) loops is indeed necessary^: 
the resummation in (|82| ) takes neatly care of this. However, the qualitative result of Ref. |2Cf| remains valid: at 
some "critical" temperature the kink solution disappears. This is not evident in the present case because of the 
dimensionality (reflecting in the temperature dependence of the propagator), nevertheless the above analysis till 



10 It should be noted that the paradigmatic example of this fact is HTL resummation which is used to cure the breakdown 
of the conventional perturbative theory for infrared momenta in thermal field theories with light bosons (especially in gauge 
theories where symmetries prevent perturbative radiative generation of masses) 
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first line of Eg. (f79|) is general and can be applied, for example, to A-0 4 in (3 + 1) dimensions: there, in the thin wall 
approximation, analytic kink solutions with spherical symmetry have essentially the same form as in (1 + 1) dimensions 
]69t . It is definitely an intriguing question if in this case the temperature at which the kink disappears coincides or 
not with the critical temperature at which the symmetry is restored. Note that both of these temperatures can be 
estimated in the present approach: for instance, the one for symmetry restoration can be read from the thermal order 
parameter in the case of homogeneous condensation (K — 0). Work in this direction is currently in progress (l9| . 



V. CONCLUSIONS AND OUTLOOK 



In this paper we have shown how to describe topological defects as inhomogeneous condensates in Quantum Field 
Theory, both at zero and at finite temperature. To this end, we have used the Closed-Time-Path formalism, allowing 
us to obtain a closed functional expansion for the Heisenberg field operator in presence of defects: the inhomogeneity 
of the vacuum condensate is introduced by means of a "shift function" f(x) which is solution of the equations for the 
physical fields (quasiparticles). 

We applied this general scheme to a specific simple model - (1+1) dimensional Xi/j 4 theory. We constructed the 
Heisenberg operator ip* in the kink sector, showing explicitly how the known classical kink solutions arise from the 
vacuum expectation value of such an operator in the classical limit. The knowledge of ip' gives us the possibility of 
calculating quantum corrections to the soliton solutions as well as Green's functions in presence of kinks. 

We have also shown how to extend this treatment to finite temperature and considered explicitly the case of high 
temperature. In this limit, we are able to calculate the order parameter in the classical limit, showing how kink 
solutions depend on temperature. This is of particular interest in the case of higher dimensionality: then the analysis 
presented here goes through in a similar way and it is possible to study the behavior of topological defects near 
the critical temperature, as well as to calculate T c . In fact, when the case of homogeneous condensation (/ = 0) is 
considered, our formulation offers a way to the study of restoration of symmetry, alternative to the traditional analysis 
based on effective potentials. 

Natural extensions of the present work includes the sine-Gordon model and the (3+1) dimensional Aip 4 theory. For 
the sine-Gordon model, which is integrable, it is of interest to investigate the correspondence between the present 
method and the Backliind transformations. We have indeed shown that in the classical limit, the introduction of the 
shift function / allows for a linearization of the equation for the order parameter. Then for integrable systems, linear 
superposition of solutions of the (linear) equation for / give rise to multisoliton solutions for the non-linear Euler 
equation: such a correspondence could indeed offer an alternative possibility for classifying soliton solutions (i.e. by 
means of the shift function /). On the other hand, we have explicitly shown that, in the case of (1+1) dimensional 
A"0 4 (which is non-integrable), it is not possible to solve in general the recurrence relation establishing the non-linear 
to linear mapping, except for a restricted class of solutions for /. Other interesting aspects of the sine-Gordon model 
for which the present approach could be useful are the finite temperature behavior |70| and the duality with the 
Thirring model both at zero |36| and finite temperature |7l| . 

For the (3+1) dimensional Xip theory with complex field - relativistic Landau-Ginzburg theory, we already found 
an interesting relation between the Green's functions for the system with vortices and the one without vortices [fj"5[ , 
allowing us to calculate a "topological" contribution to the hydrostatic pressure |ljj . Temperature-induced restoration 
of symmetry is an another important aspect which can be studied in the present framework. 

Finally, it is definitely a challenging task to extend the above formulation to the description of phase transitions in 
which defects could be created. The main problem in doing this is in the change of the set of the asymptotic fields 
when crossing different phases, which is not taken into account in the present formulation. Study in this direction is 
in progress. 
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APPENDIX A 

We show here that D[f, H — > 0] = /. To prove this we start with the definition of D[K, K— ► 0]. Let us remind that 

D[K,h^Q]=K{x)Y J J^a n b n , (87) 



n=0 



with a, b being defined in Section III. A. We shall now show that a k b l [K = /] = liffc = / = Q and otherwise (both 
k and / are positive integers). In the case when fc = Z = or Z = 0, fc^O our statement obviously holds. If k = 
but I 7^ our statement is true due to the fact that the contour integration of the function which is continuous across 
the real-time axis is zero (clearly / cannot be discontinuous as it is a solution of the homogeneous second-order 
differential equation). So let us concentrate on the remaining cases and to prove that then a k b l [K = f] =0. 

To show this let us formulate the following conjecture: 

a k b l [K}= f d 2 xd 2 yA c (x;y){s(x~y)F^[K](x)+B^[K}(x,y)} , (88) 
J c 

where both F^ k,l ^[K](x) and B^ k ' l ^[K]{x,y) are continuous functions across the real-time axis when K = f, i.e. 

Disc X0 F^[f](x) = Disc X0 B^[f](x,y) = Disc yo B^[f](x,y) = . 
To prove the conjecture (|8q ) we shall use the mathematical induction w.r.t. k. For k = 1 the conjecture is true as 

d 2 xd 2 y A c (x;y) (s(x-y)F^[K]{x)+B^[K}(x,y)\ . (89) 

From the structure of the b term is obvious that both F^' l '[f\ and B^' l '[f] are continuous across the real-time axis. 

In the next induction step we shall assume that the relation (^8f) is valid for k — n and we should prove the validity 
for k = n + 1. In the latter case we may write 

a n+1 b l [K] =a{a n b l )[K] = [ d 2 xd 2 y A c (x; y) - — 5 -— - f d 2 z x d 2 z 2 A c (z i; z 2 ) 

Jc 8K{x)5K{y) J c 

x {bW[K](z 1 ,z 2 ) + 5(z 1 -z 2 )fW[K}(z 1 )} . (90) 
Because of the assumed property of B^ n ' 1 ^ and F^ n '^ they may by written as 

BM [K](x, y) = J2 f Pq (x, y)KV{x)K q {y) P pq [b] 

P,Q 

F^[K)(x) = J2 fp(*)K p (x) P p [b] , (91) 
p 

Where both f PQ and f p are continuous across the real-time axes, -P pg [&] and P p [b] are some polynomials in b (not in 
K). Plugging ( pi] ) into ( pPf ) we obtain after some manipulations that 



a n+1 b l [K] = f d 2 xd 2 yA c (x;y){B( n+1 ' l '>[K}(x,y) + 6{x-y)F ( - n+1 ^[K}(x)}, (92) 

with an d ^(«+i,0 [/] fulfilling the required conditions. For example F^ n+1 ' l ^[K} can be found to be 

F^+^[K](x) = Ah(0) 52f P (x)p(p- l)KP- 2 (x)P p [b] 

p 

- iAn(O) [3g 2 K 2 (x)+6mgK(x)]^j^ Jj 2 Zl £>(«i)tf*Csi) . 
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During the previous derivations we have used the fact that 

d 2 yA c (x;y)r(y) = i / d 2 yG R (x,y)f p (y) = z(x) 



c 



(93) 



where z(x) is a continuous across xq. This is obvious as z{x) solves equation (d 2 + m 2 )z(x) = f p (x) with the boundary 
conditions z(0,xi) = and dz(0, x\)/dxa = (see for example [|t2| ). The second important trick which we have used 
was that 



d 2 xd 2 yA c (x,y)S(x-y)... =An(0) / d 2 x . . 
c Jc 



This closes the proof of our conjecture (|38|). As a result we may write 

a fc 6'[/] = A n (0) / d 2 xF^[f]{x) + f d 2 xd 2 y A c (x; y) [/] (x, y) = . 
Jc Jc 



(94) 



(95) 



In the previous we have applied both the fact that the contour integral over the continuous function across the 
real-time axis is zero and the identity 



A n -Ai 2 -A 



21 T ^22 



0, 



(96) 



which renders the double-contour integration zeropj. As a byproduct we obtain that only the first term in ( |87| ) 
survives, i.e. D[f, h — ► 0] = /. 

One of the main virtue of our derivation is that it immediately extends to finite temperatures because all the steps 
in the above proof may be repeated almost word by word. This mainly goes into account of equations ( [93] ) , ( pi] ) and 
(|96|) (and the arguments mentioned therein) which retain their validity even at the finite temperature level. 
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In this appendix we aim to prove the recursion relation (56). Let us first observe that for n = 2 we immediately 
obtain 



P 2 (x) = / d 2 yG R (x,y) 

J — GO 

with Pi = K. Similarly, for n = 3 we may write 
Pz{x) = i / d 2 yG R (x,y) a 



\mgPM Pi(v) + \g 2 Px{y)Px{y)Px{y) 



(97) 



\mgK 2 {y) + \g 2 K\y) 



d 2 yG R {x,y) 



|m 5 2 P 1 (y) P 2 {y) + ^3 f\(y) P 1 (y) P 2 (y) 



(98) 



During the previous derivation we have took advantage of the relation A(x, y) = — A(y, x). No assumption about the 
actual behavior of K(x) was made at this stage. As a result we may conjecture that for the general n we have 



Pn(x) 



d 2 yG R {x,y) 



mg Hy)PM + \g 2 E p M P M Pk(y) 



i-\-j—n 



i-\-j-\-k— n-f 1 



n > 2. 



(99) 



n As the identity ( p6[ ) is only based on the fact that iAn(x) = (^22(2;))* = 9(xo)iA2i (x) + 6(— xa)iAi2{x), it holds both for 
Ac(x;y) and Ac(x;y, T). Note that the latter automatically implies that 

A? a - A? 2 - A|i + Al 2 = (An + A22) 2 - (A12 + A21) 2 - 2An A 22 + 2Ai 2 A 2 i = , 



and so J c d 2 xd 2 y A(x; y)A(x; y) F(x, y) = provided Disc XQ F{x,y) — Disc V0 F(x,y) = 0. 



21 



To prove this conjecture we shall use the mathematical induction. Let us assume that Eq.(|99|) holds for n. So for 
P n+ i(x) we may directly write 



P n+1 (x) = a[P n (x)b] 



1 



d 2 yG R (x,y) I -mga 



£ Pi(y)Pj(y)b 

i-\-j=n 



1 1 



2 + off a 



n 2 T 



Piiv) Pjiv) Pk(y)b 

i+j+k—n+l 



1 f°° 

= — j d 2 yG R (x,y) 
This proves our conjecture. 



\mg ]T Pi(y)Pj(y) +\g 2 Pi(y) Pj(y) P k (y) 



i+j=n+l 



i+j+k=n+2 



(100) 



We calculate here the following integral: 



with f(z) = Ae- m ^ Zl - UZo \ We get 



APPENDIX C 



d 2 zA c (x-z)f n (z) 



(101) 



oo poo 



In = A n / dz / dzi {An (a: - z) - A 12 (x - z)} e -™»r(»i-«-o) (102) 

J ~QQ J — OO 

By performing the change of variables x— z = y and using Eq.(|52|) for An (a; — z) — A 12 (x — z) = 0(x Q — z Q )[4>(x), <f>(z)}, 
we get 



I n = - l A n e- nm ^- ux ^ dy d yi J [mJyl-yl)e nm ^- u ^ 

Z JO J-yo 



A n e - nmi ( Xl ~ux Q ) / dy e- nm "i uyo y / dw J [my y/l - w 2 ] e nm ~* v ° w 



io J-x 

where we used w = 2/1/2/0 • We then have (see |7q | formula 6.616 n.5 after analytic continuation): 



In = ~iA n - 



-iA 71 



m\/ n 2 7 2 — 1 Jo 



dy e~ nm i uyo sinh(TO2/o\/n 2 7 2 - 1) 



dy' e V.. 2 T 2 -i y ° sinh(i/ ) 



m 2 (n 2 7 2 — 1) Jo 
The last integral is tabulated (see |73], 3.541 n.l) with the result 



1^ ir(f-i)r(2) 



dxe-^sinh( X ) = -Bl^--,2) = 



4 r(f + | 



) - 1 ' 



where /i = wyu/ yjn 2 ^/ 2 — 1. The final result is then 



— -_ !_ = t A n 



?i 2 7 2 — 1 



m 2 (n 2 "f 2 — 1) /i 2 — 1 m 2 (n 2 7 2 — 1) n 2 7 2 (l — u 2 ) — 1 



= iA r ' 



> —nm'y(xi — uxo) 

m 2 (n 2 — 1) 



(103) 



(104) 



(105) 



(106) 
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So J c d 2 z Ac(x — z)f n (z) oc f n (x) provided / is an exponential solution of the Klein-Gordon equation (|34|). On the 
other hand an exponential function is the only one which has this property for all n. To see this let us reverse the 
former integral equation, i.e. 

(d 2 +m 2 )f n = (d+d-+m 2 )f n = V f n , (107) 

where r\ is a proportionality constant and d± = d/dx^ with x = (x° ± x 1 ). As a result we get for n > 1 the 
non-linear differential equation 

with £ = \J m 2 1 (n — 1) + rj/n(n — 1) which can be equivalently written as 

(0+/ + £/)(»-/ + £/) = o. 

The latter has clearly only exponential solutions. 



APPENDIX D 



In this appendix we show how the classical kink solutions emerge in the small coupling (large mass) limit. Let us put 
K — > K/g. Then, as already done in Section III.B, we may argue that since C c t is of order 0(g 2 ) it automatically 
decouples in the following reasonings. As a result we obtain: 



D[K]{x) = exp[/i<? 2 a] exp 



1 

hg< 



B 



f s 

C[K](x) = hg d 2 yA c (x,y) — — exp[% 2 a] exp 
Jc oK(y) 



hg 2 



•B 



(108) 



with B — — | J c d 2 z — ^ + ujqK 3 (z) . Keeping only the leading terms in the g — > limit (i.e. terms of order g x ) 
we obtain 



D[K](x,g^0) 



K(x) 
9 



ReSg->Q ( - exp[% 2 a] exp 



B 



hg 2 



\ _ K(x) ^ 1 , 

J 9 hw? 



B n = 



K(x) 
9 



(109) 



where the passage between second and third identity was done by means of the result in Appendix A. Analogously, 
for the C[K] term we arrive at 



C[K](x,g^0) 



9 Jc 
i 
9 



d 2 yA c {x,y) 



d yG R (x,y) 



5K(y) 
5 



Resg^o ( 9 ex-p[hg 2 a] exp 



r>9 



;B 



oc 

5K(y) Z-r „!(„ + !)! 



(110) 



Let us remark that h completely disappeared both from (]109j) and ( jL10| 
Similarly as in the case of the h limit we may write 



r ■ />oo A °° 1 

(iP f (x)) g ^ =v+ T - + - d 2 y G R (x, y) —— £ a n B n+1 
9 9 J-oo SK(y) ±^ n n\(n + l)\ 



K=f 



1 X 

V+- Rn[K}(x)\ K=f , 

9 71=1 



(111) 



where 



Ri(x) = K(x) 



Rn(x) = 



1 



[(n-2) 



112 



d 2 yG R (x 7 y)a 



n-2 



^ K 2 (y) + l -K\y) 



B 



n-2 



; n > 2. 



(112) 
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It is seen comparing Eg. (112) with Eq.(|54|) that the following recurrence relation holds: 



Rn{x) 



d 2 yG R {x,y) 



^ E Ri(y)RM + l E Ri{y)RAy)My) 



i+j—n 

The "diagonalized" recurrence relation reads 
Si{x) = K{x) 

S 2 (x)= r d 2 yG R (x,y)^- E W^O/) 



i-\-j-\-k— n+1 



n > 2. 



(113) 



i+i=2 



/oo 
d 2 yG R {x,y) 
-OO 



i-\-j—n 



i-\-j-\-k—n 



n > 3. 



(114) 



This is, as before, the functional equation of Cauchy-Marley's type. The corresponding analytical solution has the 
form S n (x) oc (Si(x)) Tl , with S\(x) = f(x) being an exponential, Fourier non-transformable and g independent 
solution of the dynamical equation (|||) . So 

S n (x) = A n f n (x) = A n e- Uon ^- x ^ . 
Incorporating this solution in the recurrence relation (|l 14j) , we are led to the conclusion that 



(n 2 - 1) I 2^o 5 ~ AlA ' 3 + 2w 2 ^ 



j4.9 -/I -j j4 



(115) 



By analogy with the calculation of Section III.B it is evident that the fundamental solution of this recurrence equation 
has the form A n — 2loq (^ZJ^j where s is a real constant. The final form for the order parameter in the small coupling 
limit may be written as 

(116) 



Recalling that v — ±u>o/g and choosing s = ±2wo e" o7a we can resum the series ( |ll6] ) with the result: 



for v positive and s negative, and 



(ip f (x)) g ^o = vth —i((xi-a)-ux ) 



(ip f (x)) g ^ = -vth -£-f((xi - a) - ux ) 



(117) 



(118) 



for both v and s negative. If s were positive, the expression ( |116[) for the order parameter becomes the non-solitonic 
one with cth[. . . ] instead of th[. . .]. An analogous analysis can be done for / = e ^oi{xi-x u) _ 
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